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PART  ONE 

THEORY  ANB  EXAMPLES 


I . INTRODUCTION 

In  a previous  report  [1],  it  was  demonstrated  that  for  linear 
and  circular  arrays  of  dipoles  the  radiation  characteristics  can  be 
controlled  by  using  a single  feed  and  reactively  loading  the  other 
dipoles.  By  varying  the  load  reactances,  the  antenna  beam  can  be 
steered. 

This  report  uses  a similar  technique  applied  to  a waveguide- 
backed  aperture  array.  Only  one  aperture  is  excited.  The  remaining 
apertures  are  reactively  loaded  by  placing  electrical  short  circuits 
in  waveguides  at  variable  distances  from  the  apertures  (see  Fig.  1) . 

An  important  advantage  of  this  design  is  a reduction  in  the  number  of 
direct-fed  elements  in  an  aperture  array. 

Coe  and  Held  [2]  considered  a similar  problem  which  involved 
the  dual  of  the  YAGI-UDA  linear  antenna  array  using  rectangular 
apertures  (reflector,  feed,  and  directors)  in  an  infinite  conducting 
ground  plane.  Their  work  was  concerned  primarily  with  radiation  in 
the  end-fire  direction. 

II . GENERAL  FORMULATION 

The  basis  for  this  formulation  can  be  found  in  a previous 
report  by  Harrington  and  Mautz  [3]  which  presents  a general  treatment 
of  aperture  problems.  We  first  consider  the  problem  of  a single  waveguide 
backed  aperture  radiating  into  a half-space  (see  Fig.  2).  The  Equivalence 
Principle  [4]  is  used  to  divide  this  problem  into  two  separate  regions 
as  follows  (see  Fig.  3).  The  aperture  is  covered  by  an  electric  conductor 
The  fields  in  the  waveguide  region  are  produced  by  the  impressed  sources 
M^,  and  the  equivalent  magnetic  current  M where 

M = A X E (1) 

AM 

over  the  aperture  region  with  the  aperture  covered  by  an  electric 


Fig.  2.  A single  waveguide-backed  aperture  radiating  Into 
half  space  bounded  by  an  electric  conductor. 
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(b)  EQUIVALENCE  FOR 
HALF-SPACE  REGION 


Fig.  3.  The  waveguide-backed  aperture  problem  divided  Into  two 
regions  with  equivalent  sources. 


conductor.  The  fields  In  the  half-space  region  are  produced  by  the 
equivalent  magnetic  current,  -Mj  with  the  aperture  covered  by  an 
electric  conductor.  The  fact  that  the  equivalent  magnetic  current  in 
the  waveguide  region  is  +M,  and  in  the  half-space  region  -M  ensures 
that  the  tangential  component  of  electric  field  Is  continuous  across 
the  aperture. 

Another  necessary  boundary  condition  is  the  continuity  of  the 
tangential  component  of  magnetic  field  across  the  aperture.  The  tan- 

V2 

gentlal  magnetic  field  over  the  aperture  on  the  waveguide  side, 

Is  equal  to 


St 


Wg 


(M) 


(2) 


where 

Is  the  tangential  magnetic  field  due  to  impressed 
sources 

is  the  tangential  magnetic  field  due  to  the 
equivalent  magnetic  source  M. 


On  the  half-space  side  of  the  aperture  we  have 


hs 


t 


(3) 


Note  that  h'^®(M),  and  are  all  computed  with  an  electric 

conductor  covering  the  aperture.  A true  solution  is  obtained  when 
of  (2)  equals  of  (3),  or 


+ H^®(M)  = - (4) 

Equation  (4)  is  the  defining  relationship  for  determining  the  equivalent 
magnetic  current  M. 

In  reality,  only  an  approximate  solution  of  equation  (4)  can  be 
obtained.  The  technique  used  in  this  report  for  solving  equation  (4) 
is  the  method  of  moments  [5].  At  this  point  in  the  formulation,  we 
will  extend  our  results  to  consider  the  multiple  aperture  case. 
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The  extension  from  the  single  waveguide  backed  aperture  problem 
to  an  array  of  waveguide -backed  apertures  with  one  waveguide  fed  is 

k+1 

straightforward.  Assume  k(odd)  apertures  with  the  center  aperture  (~^) 

driven.  An  electric  conductor  is  placed  over  each  aperture  and  the 
following  boundary  conditions  are  met  (see  Fig.  4) 


. k . . f 0 1 

(M  ) + I H 1 

4i  j=i  ^ij  (-h;  1 


i = 1,2,. ...k  (5) 

where  Is  evaluated  in  the  equivalent  waveguide  region  and 

(M^)4n  the  equivalent  half-space  region  of  the  1-th  aperture. 

M-l  is  the  equivalent  magnetic  current  M in  the  j-th  aperture  region. 
Now  define 

= u M.  (6) 

^ j J 

where  is  a complex  constant  to  be  determined  and  is  an  expan- 
sion function  to  be  specified.  Substituting  (6)  into  (5),  we  obtain 


V (u  M ) + I V (u  M ) =\ 


i = 1,2 , . . . ,k 


Next,  we  define  a symmetric  product 


<A,B>  = 


A • B dS 


aperture 


and  a testing  function  in  the  i-th  aperture  region.  Then,  taking 
the  symmetric  product  of  (7)  with  each  testing  function,  we  obtain 


/.<W  . (M  )>  + I V <W  (M  )> 

11  j=l  J '^Ij  J 


„i  ^ k+1 

-<W^,H^>  1 - — 


1 = 1,2 k 


1 
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Fig-  4.  The  waveguide-backed  aperture  array  problem  divided  Into 
multiple  regions  with  equivalent  sources. 
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Solution  of  this  set  of  linear  equations  determines  the  coefficients 
Vj  and,  therefore,  the  magnetic  current  . Once  has  been  found, 
the  field  and  field-related  parameters  can  be  computed. 

Equation  (9)  can  be  rewritten  In  matrix  notation  as  follows: 
Define  an  admittance  matrix  for  the  waveguide  regions  as 


[Y'fJ]  = [<-W  , (M  )>]  (lo; 

1 ^11  X 

and  for  the  half-space  regions  as 

= [<-w,,  (m  )>]  (ii: 

hs 

where  (M^)  is  the  tangential  magnetic  field  in  the  i-th  aperture 

region  geilerated  by  the  magnetic  current  in  the  j-th  aperture  region. 
The  minus  sign  Is  due  to  conventional  power  considerations.  Define  a 


source  vector 


= [<W^,  H^>]j 


and  a coefficient  vector 

5 - . (13) 

The  resulting  matrix  equation  which  is  equivalent  to  (9)  is 

The  physical  Interpretation  of  (14)  is  that  of  two  generalized 
admittance  networks,  and  [Y^*],  in  parallel  with  the  current 

source  (see  Fig.  5).  By  Inverting  (14),  we  obtain  the  resulting 
voltage  vector  which  is  the  vector  of  coefficients  which  determines  M 


= [Y^'S  + 


where 


0 


10 


[y"«] 


^wg 

11 


0 


0 


.Wg 

kk 


J 


,hs 

„hs 

11 

^12 

-hs 

„hs 

21 

• 

^22 

• 

,hs 

kl 

. Y 


hs 

Ik 


• Y 


hs 

kk 


An  Important  result  of  this  formulation  is  the  separation  of  the 

original  problem  into  two  regions  (waveguide  and  half-space)  whose 

W2  ll8 

characteristics  are  defined  by  the  admittance  matrices,  [Y  and  [Y  ]. 
Computation  of  [y'^®]  Involves  only  the  waveguide  region  and  [Y^®]  only 
the  half -space  region.  (One  can  see  that  a number  of  aperture  boundary 
value  problems  can  be  solved  by  this  equivalence  formulation  — see 
reference  [3].) 

Now  the  following  assumptions  are  made  (see  Fig.  6a) : 

1.  The  waveguides  and  apertures  are  one-half  wavelength  long. 

2.  The  waveguides  and  apertures  are  filled  with  a dielectric  e^. 

3.  Only  the  dominant  mode  exists  in  the  waveguides  and  apertures. 

4.  The  apertures  are  located  on  the  x-axis  with  the  long  dimen- 

sion in  the  y-directlon. 

5.  There  are  an  odd  number  of  apertures. 

6.  A single  expansion  function  represents  the  magnetic  current 
covering  each  aperture  region. 

In  addition  to  the  preceding  assumptions,  the  expansion  (M^)  and 
testing  (Wj)  functions  are  defined  as  follows: 


where 


■ '■jW  cos  ky 


Wj  ■ Mj  = Pj  (x)  cos  ky 


(15) 

(16) 


Pj(x)  - 


d -i“<x<d  +“  (see  Fig.  6b) 

V wL  J ^ J ^ 


all  other  x . 


Substituting  (15)  and  (16)  Into  (10)  and  (11)  we  obtain 


^11 

HWg 

11 

(cos 

ky)> 

(17) 

vhs  2 

’ij  ~ 

hs 

(cos 

ky)>  . 

(18) 

III.  ADMITTANCE  FORMULATION 

a)  Determination  of  y”®  - 

For  the  specific  problem  considered  in  this  report,  the  wave- 
guides are  taken  to  be  filled  with  a dielectric  (e^) , which  lowers  the 
cut-off  frequency  and  allows  the  dominant  mode  cosine  field  to  exist 
across  the  aperture. 


For  all  of  the  waveguides  except  the  externally  driven  one 
we  have 


ii 

- - j cot  1 y ( 

where  d (()i)  is  the  specific  distance  of  the  short  circuit  from 
the  aperture  in  the  i-th  waveguide  (dependent  on 
beam  steering  angle  (ji) 
k is  the  free-space  wave  number 
^c  cut-off  wave  number 

k-  - k'^  - k^ 

Ire 

is  the  relative  permittivity  of  the  dielectric  filler 
n is  the  free-space  intrinsic  Impedance 
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Equation  (19)  Is  identical  to  the  Input  admittance  of  a short-circuited 
transmission  line  where  the  characteristic  admittance  of  the  transmission 
line  Is  replaced  by  the  wave  admittance  of  the  dominant  mode,  Y^.  Physically, 
we  have  a waveguide  excited  at  one  end  by  an  external  field  (mutual  coupling) 
which  is  represented  by  an  equivalent  magnetic  current  source  and  short 
circuited  at  the  other  end  (see  Fig.  7) . 
k+1 

For  1 ■=  (externally  driven  waveguide) , we  have 


cos  ky. 


(cos  ky)> 

11 


= Y which  Is  the  wave  admittance  for  the 
o 

dominant  mode.  (20) 


b)  Determination  of  Y^^  - 

hs 

For  evaluation  of  the  off-diagonal  elements  of  [Y  ] , the 
equivalent  magnetic  current  sheet  in  each  aperture  region  Is  approxi- 
mated by  a filament  of  magnetic  current,  WP(x)  cos  ky,  at  the  center 
of  the  aperture  (see  Fig.  8(a)). 

Its 

For  evaluation  of  the  diagonal  elements  of  [Y  ] , the  concept 
of  equivalent  radius  (see  reference  [6])  is  used.  This  concept  simply 
equates  the  characteristics  of  an  aperture  antenna  to  that  of  a cylin- 
drical dipole  with  a radius  w/4  (see  Fig.  8(b)).  Even  though  this 
equivalent  radius  Is  determined  from  a quasl-statlc  principle,  the 
admittance  value  obtained  agrees  quite  closely  to  a time-harmonic 
derivation  which  solves  for  the  admittance  directly  (see  reference  7) . 


Now  both  diagonal  and  off-diagonal  terms  of  the  admittance 
lis 

matrix  [Y  ] are  In  a form  that  can  be  evaluated  using  a dipole  formu- 
Its 

latlon.  Y^j  is  proportional  to  in  [8  - equation  2-5] 


L/2 


sin  k 


(Z 


ij 


-L/2 


EijW 


sin 




kL 


z|) 


dz} 


Therefore 


WAVEGUIDE 


VARIABLE  ELECTRICAL 
SHORT  CIRCUIT 


Fig.  7.  The  reactive  load  specified  by  equation  (19) 
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Fig.  8.  Magnetic  current  source  approximations  used  In  the  evaluation 
of  diagonal  and  off-diagonal  components  of  ] • 


(21) 


[Y  ] 


[—  <-  cos  ky,  (cos  ky)>] 

wL 


2 * L * * 


In  the  above  equation,  the  first  factor  of  2 Is  due  to  Imaging  the 


magnetic  current  expansion  function  across  the  ground  plane.  The  factor 
2 

w arises  because  both  the  expansion  and  testing  magnetic  current  fila- 


ments are  proportional  tow.  Finally,  since  the  magnetic  field  due  to  a 


magnetic  current  Is  the  same  as  the  electric  field  due  to  an  electric 


current  with  u and  e Interchanged,  the  factor  -~z  la  needed  to  reciprocate 


the  factor  n Included  In  the  expression  for  E^j(z)  to  obtain  h" 


hs 


By 


■Ij 


Incorporating  the  solution  given  for  In  [8  - equation  2-7]  Into  (21), 


we  obtain 


yJ®  = -;^  (^)  [2  Cl(kd)  - CKu^)  - C1(V2)] 


- (J)  [2  Sl(kd)  - Sl(u2)  - Si(v2)] 


(22) 


where 


U2  = k (i/d^  + + L) 


V-  = k (7d^  + - L) 


Cl(x) 


00 

f cos  V 

" " J V 


dv 


Sl(x) 


sin  V 


dv 


w 


i = 3.  d 


1 j,  d Is  the  distance  between  the  apertures  1 and  j. 


while  all  the  other  apertures  are  short  circuited  (therefore,  acting  like 
a continuous  ground  plane) , a magnetic  current  M will  exist  in  the  presence 
of  an  electric  ground  plane  radiating  into  the  half-space  (z  > 0) . To 
determine  the  magnetic  field  at  a point  r , consider  the  following  analysis 
(see  Figs.  9(a)  and  (b)). 

By  placing  a magnetic  dipole  KJl  at  r and  using  reciprocity  in 

•'MU 

terms  of  this  field  an  1 the  original  field,  we  obtain  the  following  equa- 
tion 


H Ki  = 
m m 


M • H™dS 


(23) 


aperture 


,m 


where  lj[  is  the  magnetic  field  from  in  the  presence  of  a complete 
conductor  and  is  the  component  in  the  direction  of  of  the  magnetic 
field  at  r^  due  to  -M  in  the  presence  of  a complete  conductor.  The  mag- 
netic current  M in  (23)  is  generalized  to  the  set  of  magnetic  currents 
Uy  in  the  i-th  aperture  to  obtain 


H Ki  = y V,  <-M,  , H > 
m m " 1 1 


= 


i = 1,2 , . . . ,k 


?tn 


where  I is  the  transpose  of  the  measurement  vector 


r 


^1’  «”>’kxi  • 


(24) 

(25) 


Substituting  (14)  for  V in  (24)  we  obtain 

H Ki  = . 

m m 

Now  consider  the  following  assumptions: 

(1)  We  adjust  so  that  it  produces  a unit  plane  wave 
in  the  vicinity  of  the  origin,  that  is. 


(26) 


4 


1 _ -jcje  •’  m 

Kil  ” 4irr  ® 

m m 


(2)  Define  the  far-fleld  measurement  vector  as 


i“  = p®  = - 2 
j j 


-jk  • r 

M.  u • u e dS 

j *>y  «-m 


j-th 

aperture 


Substituting  (27)  and  (28)  Into  (26)  we  obtain 


H ® pntjyWg  Y^s^  1^1 

m 4tt  r 

m 


The  complex  power  transmitted  through  the  aperture  Is 


° J j £ S * 

aperture 


M • H dS  . 


aperture 


Since  this  transmitted  power  Is  only  dependent  on  the  tangential  com- 

lis 

ponent  of  H In  the  half-space  (H^  equation  (30)  becomes 

f 0 

P “ - I I M.  u • (H^^®(u  Mj)*dS 

t 1 4 i J JJ  ^ J 

^ aperture  ^ 


* „hs* 


y y V V y"® 

i 5 1 j ij 


" f hs 
V[Y  ] V 


l,j  — l«2,...,k 


The  gain  (ratio  of  radiation  Intensity  In  a given  direction  to 
the  radiation  which  would  exist  If  the  total  power,  p^,  were  radiated 
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uniformly  over  the  half-space)  associated  with  the  component  of  the 
magnetic  field  in  the  half-space  (z  > 0)  is  given  by 


2ir  r^  nlH 
m ' m' 

Real  (Pj.) 


By  substitution  of  P and  H into  (32),  we  obtain 

t tn 


G 

Sttti  Real(V[Y  ] V ) 


where 


a)  For  E-plane  patterns  - 


+w/2+x. 


p"*  = _ 2 J — 
E y wL 


Jkx  cos  4i  , ( , , 

e"^  dx  cos  ky  dy 


-w/2+x 


kw 

° ~k~  V L ® { “ kw  ) 

\ cos  Ip  ' 


b)  For  H-plane  patterns 


A/  H 

P^  = ■'2  7^  wC-sin  9)  ‘^°®^cos  ky  dy 

H y wL 

-\/4 


k Y L sin  e 


Note:  0 is  measured  from  the  positive  y axis  in  the  x=0  plane. 


V.  REPRESENTATIVE  COMPUTATIONS 

A computer  program  has  been  written  using  the  preceding  derived 
equations.  This  program  is  described  and  listed  in  Part  Two  of  this 
report.  For  this  section,  results  will  be  given  for  N=7  and  N=9  aper- 
ture arrays. 
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For  the  initial  reactive  loads,  those  which  resonated  complex  and 
real  equivalent  sources  were  tried.  In  the  4i  = 0“  case,  e..ch  required 
about  the  same  amount  of  CPU  time,  while  at  <j)  = 30°,  60°  and  90°,  the 
loads  which  resonate  complex  equivalent  sources  required  fewer  iterations 
to  achieve  a maximum.  In  every  case,  the  same  end  point  was  realized. 

Figure  10  Illustrates  the  maximum  gain  case  for  N=7  aperture  array. 
If  all  the  apertures  are  fed,  maximum  gain  can  be  obtained  since  we  have 
complete  control  over  the  complex  source  for  each  aperture.  The  disad- 
vantage of  this  case  lies  in  the  excitation.  A multiple  feed  network 
is  required  to  realize  the  various  complex  excitations  which  are  dependent 
on  the  beam  steering  angle  (j>. 

Figure  11  Illustrates  the  maximum  gain  that  can  be  obtained  using 
a single  feed  and  reactive  loads.  The  main  advantage  of  this  technique 
is  the  elimination  of  the  complex  rf  feed  networks  used  in  the  former 
case.  Reduced  controllability  of  the  gain  characteristics  is  the  penalty 
that  is  paid  by  using  this  form  of  excitation.  The  results  for  this  case 
were  computed  for  d^  = 0.27X  which  is  in  a higher  0 region  of  operation 
for  the  array. 

Figure  12  illustrates  the  same  number  of  apertures  but  in  a lower 
Q (greater  usable  bandwidth)  region  of  operation.  The  beamwldth  Is  ex- 
ceptionally high  at  <(>  = 30°  but  becomes  more  satisfactory  at  greater 
beam  steering  angles.  At  the  present  time,  this  result  would  appear  to 
restrict  the  useful  scanning  region  for  a practical  Q and  beamwldth. 

Figure  13  Illustrates  the  N=9  aperture  case.  The  gain  char- 
acteristics appear  satisfactory  except  for  the  back  lobe  direction 
((ji  = 180°)  at  low  beam  steering  angles. 

VI.  REALIZATION  OF  THE  REACTIVE  LOADS 

The  reactive  load  as  specified  by  equalton  (9)  is  short  circuit 
terminating  a transmission  line.  It  can  be  realized  using  a variety  of 
techniques.  One  of  the  simplest  realizations  is  a sliding  electrical 
conductor  In  a waveguide.  This  technique  would  be  advantageous  if  it 
were  desired  to  point  the  antenna  beam  in  a fixed  direction. 


50  25  25  50  15  10  5 5 10  15 

GAIN  GAIN 

(a)  c^  = 0“  (b)  <^=30* 


15  10  5 5 10  15  15  10  5 5 10  15 

GAIN  GAIN 

(c)  = 60“  (d)  .^=90“ 


Fig.  10.  Radiation  gain  patterns  for  a seven  element  aperture  array 

when  all  elements  are  driven  by  the  complex  equivalent  voltages 
which  yield  maximum  gain,  d = 0.27 A. 
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Fig.  12,  Radiation  gain  patterns  for  a nine  element  reactlvely 
loaded  aperture  array  when  only  the  center  aperture  Is 
externally  driven,  d.  = 0.35 A- 
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Fig.  13.  Radiation  gain  patterns  for  a nine  element  reactlvely 
loaded  aperture  array  when  only  the  center  aperture  is 
externally  driven,  d^  = 0.33 A. 


26 


If  It  Is  desired  to  actively  steer  the  antenna  beam,  the  follow- 
ing technique  Is  proposed.  A waveguide  to  mlcrostrlp  transition  along 
with  p-l-n  diodes  Imbedded  In  a mlcrostrlp  transmission  line  (see  Fig.  14) 
Is  one  possible  realization  for  the  variable  reactive  load. 

Due  to  the  Increased  availability  of  mlcrostrlp  components  using 
MIC  technology,  the  use  of  waveguide  to  mlcrostrlp  translations  has  become 
more  prevalent  In  the  microwave  Industry.  J.  Heuven  of  the  Philip's 
Research  Laboratories  [9]  has  fabricated  a waveguide  to  mlcrostrlp  tran- 
sition that  exhibits  Increased  bandwidth  due  to  the  employment  of  a 
stepped  ridge  waveguide  transformer.  He  claims  that  lower  reflectivity 
than  previous  designs  has  been  achieved.  The  waveguide  to  mlcrostrlp 
transition  shown  In  Fig.  14  Is  a modification  of  his  design. 

To  realize  variable  short  circuits  terminating  a transmission 
line,  the  use  of  p-l-n  diodes  for  the  active  shorting  elements  along 
with  mlcrostrlp  for  the  transmission  line  Is  proposed.  Under  zero  and 
reverse  bias,  a p-l-n  diode  exhibits  very  high  Impedance  at  microwave 
frequencies,  whereas  at  moderate  forward  currents  It  has  a very  low 
Impedance.  These  characteristics  permit  the  use  of  p-l-n  diodes  as  the 
active  electrical  shorting  element  In  the  proposed  design. 

The  entire  circuit  Including  waveguide  to  mlcrostrlp  transition 
could  be  housed  Inside  the  reactive  load  waveguide.  In  addition,  the 
waveguide  would  provide  adequate  rf  shielding  for  the  mlcrostrlp  circuit. 

For  the  lower  frequency  range  and  lower  dielectric  constants 
normally  used  for  mlcrostrlp  circuit  fabrication,  commercially  available 
p-l-n  diodes  could  be  used  In  their  cylindrical  shaped  package,  while 
for  the  higher  frequencies  and  dielectric  constants,  the  p-l-n  diode 
chip  would  have  to  be  used.  This  leads  to  the  feasibility  of  MIC  circuit 
fabrication  techniques  which  would  provide  a rather  compact  package  for 
the  entire  reactive  load.  M.  E.  Davis  [10]  has  clearly  demonstrated  that 
this  fabrication  technology  Is  at  hand. 
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VII.  DISCUSSION  AND  CONCLUSIONS 

The  computer  program  has  Its  basic  limitations  as  does  any 

mathematical  model  of  a physical  system.  One  limitation  Involves  the 

spacing  between  apertures  and  the  aperture  width  (some  kind  of  falling 

as  In  a thln-wlre  dipole  array  problem).  This  limitation  Is  most  prevalent 

In  the  high  Q operating  region  of  the  aperture  array.  Numerically,  the 
HS 

matrix  Real  (Y  ) used  In  the  gain  calculations  Is  numerically  not  posi- 
tive definite  (see  Table  1)  and,  therefore,  the  gain  becomes  unbounded. 

The  high  Q region  of  operation  is  usually  not  the  most  useful  operating 
point.  However,  a lower  Q region  may  yield  an  undesirable  beamwldth  at 
the  lower  beam  steering  angles.  Therefore,  an  operating  point  in  between 
appears  to  be  the  most  promising. 

While  the  gain  is  not  a function  of  aperture  width,  the  aperture 
admittance  Is  related  to  it  (see  equation  22).  Use  of  the  thin  dipole 
approximation  for  calculating  the  aperture  admittances  restricted  us 
to  using  very  narrow  apertures.  These  small  widths  yielded  aperture 
admittances  and  ultimately  reactive  load  values  (which  are  of  the  same 
order  of  magnitude)  considerably  less  than  the  dominant  mode  wave  ad- 
mittance, Y^.  Tables  2 and  3 show  the  reactive  load  values  required  to 
obtain  the  data  for  the  cases  Illustrated  in  Figs.  11  and  12.  Y^  for 
these  two  cases  Is  1.97  mmho,  which  corresponds  to  a dielectric  con- 
stant of  “ 1.55.  Therefore,  due  to  the  large  difference  in  magnitudes 
between  Y^  and  the  aperture  admittances,  a mismatch  occurs  at  the  wave- 
guide aperture  to  half-space  botmdary.  In  addition,  Tables  4 and  5 depict 
the  short  circuit  distances  required  to  realize  the  reactive  loads  given 
in  Tables  2 and  3 for  the  (|)  = 0°  direction.  The  small  distances  required 
would  make  the  physical  realization  suggested  in  the  previous  section 
difficult. 

The  analysis  so  far  has  considered  only  the  lossless  case.  In 
reality,  loss  would  be  present  in  any  physical  realization.  If  loss  were 
added,  the  thin  dipole  approximation  would  yield  a positive  definite 
half-space  admittance  matrix  which  would,  in  turn,  allow  greater  aperture 
widths.  Table  6 shows  the  N=7,  d^  - 0.27X  case  considered  in  Table  4 but 
with  1%  loss  added.  This  allows  a more  reasonable  aperture  width  of  O.IX 
Figure  15  Illustrates  the  gain  patterns  in  this  lossy  case  (compare  Fig.  11 
with  Fig.  15  to  note  the  reduction  in  gain) . 
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Table  1. 


Eigenvalues  for  Real 


(Y^®)  — N=7, 


0.25)i 


APERTURE  WIDTH 


EIGENVALUE 

No. 

W = O.OIX 

W = 0.02X 

W = 0.03X 

W = 0.04X 

1 

0.86  X 10”^ 

0.17  X 10"^ 

0.26  X 10“^ 

0.34  X 10"^ 

2 

0.80  X 10“^ 

0.16  X 10“^ 

0.24  X lO"^ 

0.32  X 10"^ 

3 

0.63  X 10"^ 

0.13  X 10"^ 

0.19  X 10"^ 

0.25  X 10"^ 

4 

0.49  X 10“^ 

0.98  X lO"^ 

0.15  X 10”^ 

0.20  X 10"^ 

5 

0.94  X 10"^ 

0.19  X 10"^ 

0.28  X 10"^ 

0.38  X 10“^ 

6 

0.49  X 10“^ 

0.98  X lO"^ 

0.14  X 10~^ 

0,19  X 10“^ 

7 

0.63  X lO"® 

0.22  X lO"^ 

-0.31  X lO”^ 

-0.10  X lO"^ 

Table  2.  Reactive 

load  values 

(mmhos)  for  a i 

seven  element 

aperture 

array  for 

= 0.27X  and  W 

= O.OIX  (see 

Fig.  11) 

APERTURE 

NO. 

<|)  = 0** 

(fi  = 30° 

<t>  = 60° 

<j>  = 90° 

1 

-0.007  my 

-0.008  my 

-0.011  my 

-0.006  my 

2 

-0.026 

-0.029 

-0.033 

-0.033 

3 

-0.034 

-0.037 

-0.032 

-0.032 

5 

-0.036 

-0.032 

-0.037 

-0.032 

6 

-0.026 

-0.030 

-0.026 

-0.033 

7 

-0.057 

-0.027 

-0.021 

-0.006 

Table  3.  Reactive 

load  values 

(mmhos)  for  a i 

seven  element 

aperture 
Fig.  12) 

array  for  dj 

• 

= 0.35X  and  W 

= O.lOX  (see 

APERTURE 

NO. 

<|i  - 0“ 

<t>  = 30° 

<(.  - 60° 

(J.  = 90° 

1 

-0.121  my 

-26.4  my 

-0.235  mU 

-0.083  mU 

2 

-0.942 

-0.108 

-1.36 

-0.635 

3 

-0.442 

-0.363 

-0.191 

-0.345 

5 

-0.507 

-0.689 

-0.858 

-0.345 

6 

-0.466 

-0.842 

-0.227 

-0.635 

7 

-0.284 

-11.3 

-0.333 

-0.083 
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Fig.  15.  Radiation  gain  patterns  for  a seven  element  reactively 
loaded  aperture  array  when  only  the  center  element  is 
externally  driven,  d^  = 0.27  A,  and  a loss  of  1%  is  added. 
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Table  4.  Short  circuit  distances  for  reactive  load  realization 

for  ^ = 0°,  dj  = 0.27X,  W = O.OIX,  and  freq.  = 2.9  GHz. 


APERTURE 

NO. 

Load 

(e^  = 1.1) 

d^(<ti)  (e^  = 

1.55) 

1 

-0.007  mU 

81.4  nnn 

34.8  mm 

2 

-0.026 

80.2 

34.6 

3 

-0.034 

79.7 

34.5 

5 

-0.036 

79.6 

34.5 

6 

-0.026 

80.2 

34.6 

7 

-0.057 

78.3 

34.3 

Table  5. 

Short  circuit 

(()  = 0®,  dj  = 0 

distances  for  reactive  load 
.35X,  W = O.lOA,  and  freq. 

realization 
= 2.9  GHz. 

for 

APERTURE 

NO. 

^Load 

d^(<>)  (e^  = 1.1) 

d^(<!>)  (e^  = 

1.55) 

1 

-0.121  m{5 

74.5  mm 

33.5  mm 

2 

-0.342 

61.7 

31.1 

3 

-0.442 

56.6 

30.0 

5 

-0.507 

53.5 

29.3 

6 

-0.466 

55.4 

29.7 

7 

-0.284 

64.8 

31.7 

Table  6.  Short  circuit  distances  for  reactive  load  realization  for 
4,  = 0®,  dj  = 0.27A,  W = O.lOA,  loss  added  = 1%,  and 

freq.  = 2.9  GHz. 


APERTURE 

NO. 

®Load 

d^(<|))  (e^  = 1.1) 

d^(<),)  (e^  = 1.55) 

1 

-0.075  mU 

77.2  mm 

34.0  mm 

2 

-0.435 

56.9 

30.0 

3 

-2.622 

16.1 

14.3 

5 

-0.380 

59.7 

30.6 

6 

-0.372 

60.1 

30.7 

7 

-0.208 

69.2 

32.5 
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In  this  report,  a fabrication  technique  Is  suggested  for  the 
physical  realization  of  the  reactive  loads.  In  some  respects,  this 
reactive  load  realization  might  be  considered  more  complex  than  a 
variable  ferrite  phase  shifter  and  on  equal  par  with  switched  delay 
line  and  p-l-n  reactive  diode  phase  shifters.  However,  reduction  in 
the  number  of  direct-fed  elements  in  the  reactive  load  case  certainly 
is  advantageous. 

Work  should  be  continued  to  minimize  the  180°  back  lobe  at  the 
lower  beam  steering  angles  for  higher  N and  higher  Q-mode  aperture 
arrays.  Also,  some  attention  should  be  given  to  a possible  composite 
array  using  phase  shifters  and  reactive  loads.  Phase  shifters  could 
be  used  for  the  major  beam  steering  and  variable  reactive  loads  used 
to  obtain  other  desirable  properties  such  as  minimizing  unwanted 
sldelobes. 
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PART  TWO 

COMPUTER  PROGRAMS 

I . INTRODUCTION 

A computer  program  used  to  generate  the  figures  In  this  report 
along  with  a numerical  example  is  presented  in  this  part  of  the  re- 
port. This  program  consists  of  the  subroutines:  YHSP,  SICI,  PATHSP, 

MAXGCV,  MAXGRV,  LINER,  BLOADC,  BLOADR,  FUNCTA,  FUNCTB,  LINEQ,  EIGEN, 
and  the  calling  MAIN  program.  Each  subroutine  computes  one  function 
only  and  can  be  used  or  not  used  if  that  specific  output  is  or  is  not 
required. 

II.  ADMITTANCE  MATRIX 

The  subroutine  YHSP  (N,  X,  W,  L,  YHS)  computes  and  stores 
columnwise  in  YHS,  the  elements  of  the  admittance  matrix. 

Ci(kd)  -Ci  (u  ) - Ci(v  )] 

13  r|iT  L i.  / 

- (f>  [2  Si(kd)  - Si(u„)  - Si(v.)]  (1) 

riTT  L t i 

where 

U2  = k ( i/d^  + + L) 

V2  = k ( yd^  + L^  - L)  . 


There  are  N apertures.  X(I)  is  equal  to  kX^  where  k is  the  propagation 
constant  and  is  the  coordinate  of  the  I-th  aperture,  w is  the  aperture 
length.  The  IJ-th  element  of  the  port  admittance  matrix  is  computed 
inside  nested  DO  loops  10  and  11  and  is  stored  in  both  YHS (I  + (J-1)  * N) 
and  YHS(J  + (I-l)  * N) . 
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Minimum  allocations  are  given  by 

COMPLEX  YHS  (N*N) 

DIMENSION  X(N) 

The  subroutine  SICI (SI,CI,X)  stores  the  sine  and  cosine  integrals 


Cl(x) 


OO 

r 


cos  V 

V 


dv 


X 


f sinv, 

Si(xj  = — - — dv 

0 


(2) 

(3) 


In  SI  and  Cl  respective.  SICI  Is  the  subroutine  SICI  in  the  IBM 

* 7T 

Scientific  Subroutine  Package  with  SI  replaced  by  SI  + y. 


IBM/360  Scientific  Subroutine  Package  (360A-CM-03x)  Version  III, 
Programmer's  Manual,  page  370. 
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SUBf- JLITJ  -J£  ynSP(N,X,W,L,YHS) 

CjMPLcX  YhS(Sy9)tU 
DIMENSIC  N X(9) 

PL  AL  L,L  1 
PI=3. 141 593 
cT A=376. 730 
L = L/^!. 

L 1 = L='2.«‘  P I 

T = 2. 

U=<C.,  1.  » 

L2=L1’>L1 
Tl  =T  <'L  1 
TL2=TLoTL 
Ci=rcs(L  1 ) 

SN  = S1N(L  1 ) 

Cl=w/(  Pi  •cTA’‘L<'Sf.«Sr'l) 

C2=SN*CS 

C3=.5*CLS(TL) 

JZ=^ 

DO  10  1 ,N 

J1  =J 

00  11  I = 1 , J 
J3=J2+I 

IF(I-J)  12,13,12 

13  IFM-i)  14,15,14 

14  YHS( Jl)=  YHS( 1 ) 

GC  TG  11 

15  D=(W^'2.=  PD/4. 

02=C»D 

GO  'TO  16 
12  X j = x( n- <( jj 
02=XJ»XJ 
0=SCkT(D2  ) 

16  Sl=SOkT< 02+L2) 

S2=SQPT( 02  + TL2  ) 

V1=S1+L1 

ai=D2/Vl 

U2  = S2+-tl 
V2=C2/U2 

CALL  5lC  f (SL),CD,U» 

CALL  SIC  I (SOI  ,C'Ji  ,U1  ) 

CALL  SIC  I (SV1,CV1,V1) 

CALL  SIC  I (SU2,CU2,U2I 
CALL  SIC  I (SV2,CV2,V2) 

YHS1  = C1-  (C<;*  ( SU2-SV2-T«(  SV  1-SUl ) ) -C  3 ( 2 U 1 -CD+LV  1 ) - lU2-L  ) -CM  1 + 
1T»-CC-CVI  ) 

YHS2=C  l'^  ( C2*  ( CV2-C02«-T*-(  r,  Vl-CU  1)  ) -C  3 ( T*»  I S U 1 -SU+ 5V  1 ) - S S V«;  ) - j J 1 + 
1T^«SD-SV1  ) 

YHS(  J1  J=  YHS1-U’*YHS2 
YHS(J3)= YHS(Jl) 

J1=J1+-N 
11  CONTINUf 
J2=J24-N 
1.  CJNTINUL 
PE  TURN 
END 


SUBHOUTl  NE  SI  cl  t SI  ,f;  1 ,/) 

Z=ARS(X) 

IF(Z-A.) i ,1,^ 

I Y=  (^.-Z)  «■- (4.  + Z » 

3 SI=X-»(((((l.75  3l*lZ-'»®Y4^l.  56d  VcJ8E  - 7 ) -‘Y<- 1 . 3 loot  -0  ) » Y +0  . v 3 SonC  - 
1*Y«- 1.  v6^SB^E-2  ) *-Y4-A.  1 ) 

C I = ( ( b.7  7215>6L-1  + ALv.M  Z))/Z-Z*^(n((l.3  86Si35i--lc<-Y+l.Do‘t‘i->oL-o)’V 
1 + i .7ZS7S^c-6  ) oY*  I . l^S'^Si'yF-'t ) *Y+^t  . V9v.  V^:t-3)*Y«-l.ji3j>;«E-l)  )'Z 
(iE  TLRN 
A S I = S I N ( 7 ) 

Y = CCS(  Z) 

Z=A./Z 

'J=(<((((((A.046OOVc-3»Z-2.279  1A3£-2)~Z<-S.3  13U/Ct-^J'fZ-7.^nlf‘t2.-- 
l*Z-«-4.9ci7  7 16E-2  )*^Z-3.  3:><i51  :<E-3)*=Z-^;.314617E  - ^ t «Z-  i..i:>4  9^>3L-5)''Z 
2*-0.2S00l  lt-2)  "■Z+?.5j39d9E-lJ 
y = (((((({((  - 5 . 1''  H699E-3>"  Z + 2.  d 19  179L-2)  *Z-5 . 5 :>7^03L- t ) *^Z 
U7.90203  4r-2)*Z-H.9cU4  1bE-2)*Z-7.'>4S35oi-3)--Z+Z.6jld9  3E-2)'/- 
2-3. 764'-,  Jt-4)  *Z-Ji.  12  24  13E-2)  »Z -0.646  4411  - 7 ) Z+ ^ . 6 c o ..  C-  . t - i 
Cl  =Z*(SI  ’‘V-Yx'U) 

S I =-Z<-  (£  I *IJ+Y*'V  ) +1.670796 

RETURN 
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The  output  is  stored  in  PME(N,J)  and  PMH(N,J)  for  <))  = (J  - 1)*  ^ , 

Jo 

J = 1,2,  ... , 37. 

The  input  variables  are  the  same  as  YHS. 

Minimum  allocations  are  given  by 

COMPLEX  PME(N,37),  PMH(N,37) 

DIMENSION  X(N) 


I 


1 


SUHROLTU!-  TMSR( 'nI  , X ,1  , W f MMP  , PMM  ) 
L ex  ( > t 3 M , {'y  , 37  J , S,  SI , U 

iJIMeNSIlJiM  X ( ) ) 

RRAL  L,lI 
P(  = 3.  '^l7S;3 
u = ( n . , I . ) 

L 1 = L«<;.X'P  I 
1 = W*  PI 
DEI =P  1/36 

C 1 ( P . vw/L  ) 

Cc'  = - 1 .^Cl 
DO  10  I = 1 , N 
DO  11  J= I , )7 
PH^  ( J-1  )^l)h  I 
LS  = Cr  SI  PH) 

SN=S  IMPH) 

S=  00  S ( X ( I ) - t :>  ) ♦ S 1 M X I I)  C b ) 

1 M I . Eg.  19)  P,)^-  I 1 , J ) -t  IX^S 

1 F ( I .(-Q.  19 ) ;o  nj  1 

PMf  ( I , J ) =C  IM  K i^Cb  ) / (Wl^OS  ) 

12  COM  INOE 

IF  ( I .1.7.1  ) P3H(  !,.))- PMH  I 1,J) 

I F ( I . 0 7 . 1 ) (,  I r 1 11 

IF(J.FQ.i)  P3h( I ,J) 

If  (J.LU.3  7)  PMiK  I ,J)=Cl 
I F ( ( J.FO.  1)  . )R.  U.FU.  37  n i.O  T!'  11 

Sl=CnS(PI  <=Cb/2.  f/s-'J 
PMHI  I , J) =C3-S I 
11  CoMIMUt 
10  CONUNUL 
PE  ruPN 
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IV.  MAXIMUM  GAIN,  COMPLEX  AND  REAL  EQUIVALENT  VOLTAGE  SOURCES,  ALL 
APERTURE  FED 


(The  derivation  of  the  formulas  found  in  this  section  can  be 
found  in  reference  [11].) 

The  subroutine  MAXGCV(N,PMER,YHS ,VCMR)  computes  the  maximum 
gain,  GC,  using  complex  equivalent  voltage  sources  for  excitation. 
These  sources  are  outputted,  VCMR,  so  that  BLOADC  can  compute  the 
reactive  loads  needed  to  resonate  these  sources  for  a starting  point 
for  the  optimization  program. 


[Ghs^-1 


(6) 


VCMR  = [G^®]“^  P™*  (7) 

1 

In  the  program,  the  constant  — — has  been  changed  to  - — because 

oirr)  oirri 

the  measurement  vectors,  p'",  have  been  multiplied  by  k to  make 
them  insensitive  to  absolute  length. 


The  subroutine  MAXGRV(N,PMER,YHS,VRMR)  computes  the  maximum 
real  gain,  GR,  restricting  the  array  excitation  to  real  or  equiphasal 
equivalent  voltage  sources.  These  sources  are  outputted,  VRMR,  so 
that  BLOADR  can  compute  the  reactive  loads  which  resonate  these  sources. 


2 

GR  = -I—  Real(p”)  [G^®]“^  {Real(P™)  + c Imag  (P™)}  (8) 

oTTf) 

VRMR  = [G^®]"^  [Real(?")  + c Imag  (P™)  ] (9) 


where 


c = - a + 


(a^  + 


1) 


1/2 


Real(P^)  [G^®]  ^ Real(p”^)  - Imag  (p”*)  [G^®]~^  Imag  (p"*) 
2 ReaKP”)  [G^®]“^  Imag  (?“) 


40 


[ 

t 

i 

f 


Note:  There  are  two  possible  real  solutions  for  c which  make  GR  stationary. 

The  maximum  gain  is  obtained  by  using  the  larger  of  the  two  values. 

Minimum  allocations  are  given  by 

COMPLEX  PMER(N),  YHS(N*N) 

DIMENSION  G(N*N),  VI (N) ,V2 (N) ,A1(N) ,A2 (N) ,VRMR(N) 

The  subroutine  LINER(LL,C)  is  an  inversion  subroutine  for  real  matrices. 
Minimum  allocations  are  given  by 

DIMENSION  LR(N),C(N*N) 


i 

I 

I 


1 


I 
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SUBPOUTl NE  MAXGC V( N, PMfcR , YHS , VCMh ) 

COMPLtX  CCNJGf GC 1 (9) ,PMEP (9) , VC«P (9) ,YHS<99) 

DI  MENS  IC  N 0(99) 

PI=i. 141 993 
£TA=376.  730 
NN*N*N 

C1=1./(«.*PI*ETA) 

00  1 J=1 »NN 

G( J) =REAL (YHS( J) ) 

1  CONTINUE 

CALL  LINEP(N,G) 

DO  3 K = 1 ,N 
GC 1(K)=( 0.,0. ) 

VCMP(K)= (0.,0.  ) 

3 CONTINUE 
Ji  = l 

DO  4 K=1,N 
DO  5 J=1,N 

GCKK  )=GC1(K)  + PME«(J)*‘G(  Jl) 

J1=J1+1 

5 CONTINUE 

4 CONTINUE 
GC=C. 

DO  6 K=1 ,N 

GC  = GC^REAL(GC1(K)*‘CjNJG(PMEP  (K))  ) 

6 CONTINUE 
GC=GC*C1 
Jl  = l 

DO  7 K = 1 ,N 
DO  8 J=1 ,N 

VCMR(K)  = VCMR(K)+G(J1  I-S'CJNJGI  PMEK(  J)  ) 

J1  = J1+N 
8 CONTINUE 
J1=K*1 

7 CGNT INUE 
WPITE(3,10)  GC 
WRlTt(3fll)  I VCMH( J) ,J=1,N) 

10  FQPMAT(/3X, 'COMPLEX  GAIN  =',1F5.2) 

11  FORMAT (/ IX,  • complex  EJUI  VALENT  VuLTAGE  F UP  MAXIMUM  COMPLEX  GAIN* 
1// (3X,5E 14.7) ) 

RETURN 

END 

subroutine  MAXGPV(N,PM£.<,YHS,VkMF:) 

COMPLEX  PMER(9),YHS(99) 

DIMENS luN  G(99) , VI ( 9) , VZ(9) , A 1(9) ,A2 (9) ,VRMP (9) 

PI=3. 141593 
ET  A=376.  730 
NN=N*N 

C1=1./(8.*PI*ETA) 

DO  1 J=1 ,NN 

G(  J) =PEAL (YHS ( J) ) 

1 CONTINUE 

CALL  LINER(N,G) 

DC  2 1=1  ,N 

Vl(  I)=REAL(PMER(I  ) ) 

V2(  I )=AI MAGI PMER  ( I )) 

2 CONTINUE 
DC  3 K=1  ,N 
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VkM^(k  )=  j. 

(K ) =J. 

A?  CK  ) =^0. 

J3  = (K-1) 

Du  4 1=1  ,N 
Jl=J3*l 

A1  (K)=A1  ( Kl«-u(  J I )-‘Vl  ( 1 ) 

A2 (K ) =Ac  ( M IJi  1 >V2 ( I > 

4 CtiNT  INUe 
3 CGNTINUl 

51  = C. 

52  = C. 

S3=C. 

DO  5 K=1  f h 
S1=SH-V1  (K)«A1  {<) 

S2  = S2  + V^  (K)  vA2  (X  ) 
b3  = S3«-Vl  (K)*A2  (K  » 
b CoNTINUl 
S4=Sl-S2 

IF(S3.£g.a.»  S3=S4='l  .F-14 
A=S4/ (2. ♦S3) 

SA=SQRT(  A*A+l.  ) 

C=-A^SIG,’'nSA,S3) 

GR=r  I*  »S  UC*S3) 

WRITF(3,20)  GP 

23  FORRAH/ 3X, 'MAXI '3UM  RtAt  GAIN  =',1F5’.2) 

J3  G 1=1,  N 

vr  WR  ( n = Ain  ) -xc*  A2(  I ) 

9 CGNTINUc 

WK  I1F(  3,  2 1)  ( VFNR  ( I)  ,I  =1  ,N) 

FORMAT  </ JX, 'REAL  EQUIVALENT  VCL'^AGt  F^K  MAXIMUM  kLAL  GA  I N • / / ( .- X , 4 
114.7)) 

KcTURN 

END 

SOPRCUTI NE  LINCP (LL,C) 
dimension  LP(3r.)  ,C(99) 

DU  ?)  1=1, LL 
LR( 1 )=I 
2j  CGNTINUL 
Ml  =C 

DO  Iti  M=l,LL 
K = M 

DJ  2 I=M,LL 

K1=M-H 

K2=M1+K 

IF ( ABSIC ( K1 ) ) -AGS <C(K2 ) ) ) 2,2,6 
6 K=  1 

2 CCNTINUL 
LS=LR(M) 

Lk<M)=LF  (K) 

LR (K )=LS 
K2=mi,.K 
STCF=C(K  2 ) 

J1  =c 

DU  7 J=1,LL 
K1=JWK 
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K2*JUM 

sro=c (Ki ) 

C(KU=C(K2) 

C( K2)^STQ/ST0R 
J1=JH-LL 
7 CONTINUE 

C(K1)=1./STCR 
OG  li  1=1, LL 
IF(I-M)  12,11,12 

12  K1=M1*-I 
ST=C<KU 
C( Kl) =0. 

J1=C 

DC  IJ  J=1,LL 

K1=J14-I 

K2=J1+M 

C(K1)=C(K1»-C(K2 
J1=J1  + LI 

10  CONTINUE 

11  continue 

M1=K1+LL 
18  CONTINUE 
J1=C 

DO  <5  J = 1,LL 
IF(J-LRIJ))  14,8,14 
14  LR J=LR (J I 

J2=(LRJ-1  )’>=LL 
21  DC  13  1 = 1, LL 
K2=J2-H 
Kl=Jl+» 

s=cn 

C(K2I=C( Kl) 

C(K1)=S 

13  CONTINUE 

LP (J)=LP (LRJ) 

L?  (LRJ)=LRJ 
IFU-LR(J))  14,8,14 
a J1=J1+LL 
9 CONTINUE 
RETURN 
FNO 


J 
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V.  RESONANT  LOADS 

The  subroutine  BLOADC(N,RFP,BLOAD,VCMR,YHS)  stores  in  BLOAD, 
the  reactive  loads 

N 

Imag  [-  I Y VJ  - ReaKY^^  *Imag(V  ) 

BLOAD(I)  = ^ (10) 

Real  (V^) 

for  I = 1,...,N 

which  resonate  the  complex  equivalent  voltages  VCMR.  NFP  is  the  aperture 
feed  port  number. 

Minimum  allocations  are  given  by 

COMPLEX  VCMR(N),  YHS(N*N) 

DIMENSION  BLOAD (N) 

The  subroutine  BL0ADR(N, NET, BLOAD, VRMR.YHS)  stores  in  BLOAD, 
the  reactive  loads 

BLOAD(I)  = ([b’"®]  * (11) 

Vi  O 

(where  ( [B  ] * V)  denotes  the  I-th 

^ hs  ->■ 

component  of  the  column  matrix  [B  ] * V) 

which  resonates  the  real  equivalent  voltages,  VRMR. 

(Derivation  of  this  formula  can  be  found  in  reference  [11].) 

(Note  that  it  can  be  easily  shown  that  the  gain  does  not  depend  on 
V+1 

BLOAD(I)  where  I = — ^ (driven  aperture).) 

Minimum  allocations  are  given  by 

COMPLEX  YHS(N*N) 

DIMENSION  BLOAD (N) , VRMR (N) 


f 
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r 


I 

, SUBPOUTInE  BL0ADC(N,NFP,BL0AD,VCMR ,YhS) 

I CO»^PLEX  Ul,  VCMR(9)  ,YHS  (S^9> 

OIHENSILN  BL0AD(9| 

Jl=l 

DO  1 l=ltN 
J2  = I 

Ul  = (0.  .0 . I 
DO  2 K=1 ,N 
IF (K.EQ. I ) GO  TO  3 
Ul=tl-YHS(J2 I»VCMk(K) 

3 CONTINUE 
J2=J2*N 

2 CONTINUE 

BLOAOI  I)  = (AIMAGIUl  )-RE  AL(  YHS  (Jl)  I^^AI  MAGIVCMP  ( I ) ) ) / RE  AL(  VCMP  ( I ) » -/  I 
IMAGIYhSI Jl) ) 

J1=J14N+1 

1 CONTINUE 
BL0AD(NFP)=0. 

WRITE(3,4)  (BLOAOm  ,I=1,N) 

4 F0RMAT(/3X,' BLCAD-PESONATI NG  COMPLEX  SGUPC E S * / ( 3 X , 5E 1 4. 7 ) ) 

RETURN 
END 

SUBROUT  I NE  BLGADK (N , NFP t BLOA D t VF  MR , YhS ) 

COMPLEX  YHSI99) 

DIMENSION  BLUADI 9) ,VRMR(  9) 

DO  I 1=1, N 
J1  = I 
U1=C. 

DC  2 K=1,N 

U1=LU  (A  IMAGIYHSI  Jl)  ))*VPMRIK) 

Jl=Jl+N 

2 CONTINUE 

BLOAGI I) =-Ul/VRMP  I I ) 

1 continue 

BLOADINFP )=0, 

WRITEI3,3)  IBLOADI I ) ,l =l ,N) 

! 3 FOBMATI/ 3X,‘BL0A0-RESCNATING  REAL  SO UP Ct S * / I 3X , 5t 1 4.  7 ) ) 

RETURN 

' END 

I 

i 

t 

I 

I 

f 

I 

) 

I 
I 

1 

1 

I 
1 
I 
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VI.  GAIN  PATTERNS 

The  subroutine  FUNCTA(N,GAIN,PMER,V,YHS)  computes  the  gain  (all 
apertures  excited)  using  the  following  equation 

GAIN  = r (12) 

STrn  Real(V[Y^®]  * ^*) 

“V 

for  a given  measurement  vector  PMER  and  excitation  vector,  V. 

Minimum  allocations  are  given  by 

COMPLEX  PMER(N) ,V(N) ,V1(N) ,YHS(N*N) . 

The  subroutine  FUNCTB(N,NFP,BLOAD,GAIN,PMER,YHS)  computes  the  gain 
using  formula  (12)  with  the  following  modification.  The  excitation  is 
that  of  a unit  current  source  for  a single  aperture  (NFP)  and  the  other 
apertures  reactlvely  loaded  (BLOAD) . Therefore, 

V = [Y^®  + .1  BLOAD]  * I (NFP)  (13) 

where 

I (NFP)  = 


Minimum  allocations  are  given  by 

COMPLEX  PMER(N) ,V(N) ,V1(N) ,Y(N*N) ,YHS (N*N) 

DIMENSION  BLOAD (N) ,G(N*N)  . 

The  subroutine  LINEQ(LL,C)  is  used  in  FUNCTB  to  invert  a complex 
matrix. 

Minimum  allocations  are  given  by 

COMPLEX  C(N*N) 

DIMENSION  LR(N  ) 


0 

1 

0 


SUBPOUT  I Nc  PUNCT^.  ( N , GA I \ , PMt  P , V , Yh  S ) 

COMPLEX  CCNJG,  PMEP  (V),  V(  9 ) ,V1  YHS(  ) 

£TA=37o. 730 
PI =3. 141 S93 
CI=1./(8.*PIoETA) 

V2=(  J.,G  . » 

DO  1 J=1 ,N 
VI ( J) =(G  . ,0.  ) 

J3=( J-l) •N 
DO  2 I =1  ,N 
J4  = J3«-I 

Vl(JI=Vl  ( J)-*-CUNJG  (Vi^S(  J4>  »*CCNJG(V(  I M 

2 continue 

V2=V2*Vl  UI^V  ( J> 

I  CONTINUE 
DEN  = REAL  5V2) 

V2=(  0,  ,C  . ) 

00  3 J = 1 ,N 
V2=V2+V( J»»PMEP( J) 

3 CONTINUE 

GAIN=(Cl*V2fCCNJG(  V2U  /DEN 

RETURN 

END 


SUBROUTINE  FUNCTb  (N.NF  P,  B i.  UA  D , GA  it, , P M rK  , YH  S ) 
COMPLEX  CCNJG, PMEP (0) , ul , V ( 9 > , V 1 ( 9 J , V2 , Y ( 9 9 > 


UIMENSICN  BL0AD(9),0(99) 
ETA=  376.  730 
PI =3. 141 593 


Cl=l./(6.*PI'*ETA) 

Ul  = (0,  ,1  . ) 

NN=M*N 
Nl=N+l 
DO  I J=l,NN 
Y( J)=YHS I J» 

J) =PEAL (YHSI J) ) 

1 continue 

J3  = -N 

DO  2 J = 1 ,N 
J 3 = J3  + NI 

y ( J3) =Y( J3)+U1*HL0A0( J ) 

2 Continue 

CALL  LINEGIN.YJ 
I 3= (NFP-  I )»N 
00  3 J=l,N 
14=1 3+J 
V( J»=Y(I4> 

3 CONTINUE 
V2=(G.,i'  . ) 

00  4 J=l,N 

VI (JI  = (j  . ,u.) 

J3=( J-l) 

DO  5 1=1  ,N 
J4=J3*I 

VM  J)=Vl  ( JM-G(  J4)*V(  I) 

5 CONTINUE 

V2  = V2^V1  ( J)*CuNJG(V( J) ) 


,YhSl99) 


! 


I 

1 


CJNT  INUt 

DEN  = REAL(V2)  ^8 

V2=(0.,J.  ) 

Da  7 J = l ,N 

V2  = V2-frV(  J)*PMER(  J) 

7 CONTINUE 

GAIN=<C1»REAH  V2«'COUJG(V2)  ))  /DEN 

RETORN 

END 

DIHENSICN  LR(5d» 

DO  20  1=1, LL  i 

LRII)=I  1 

20  CONTINUE 
Mi=r 

DO  18  M=1,LL 
K = K 

<2=M1+K 

Sl=AfiS(REAL(C(K2)  n + A&S(AIMAG(C  (K2))  » 

DO  2 I=M,LL 
K1=M1-H 

S2=ABS(E  EAL (C ( K1 ) ) ) ♦ AB S ( A I MAG ( C ( K 1 1 ) ) 

IF(£2-S1  ) 2,2 ,6 

6 K=  I 
Sl  = £2 

2 CONT INUE 
LS=LR(M> 

LR (R l=LR (Kl 
L0(K)=LS 
K2=M1+K 
STCR=C(K2  ) 

J1=C 

DO  7 J=1 ,LL 
K1=J1<-K 
K2=JH-M 
STC  = C(K1  ) 

C(K1)=C(K2> 

C(K2)=ST0/ST0P 

Ji=Jl+LL 

7 CGNT  INUE 

K1=M+M  i 

C( Kl)=l. /STOR  I 

DO  11  1 = 1, LL  I 

IF(I-M)  12,11,12  1 

12  K1=M1+I 
ST=C(Kl) 

C(K1)=0. 

J1=0  j 

DO  10  J=1,LL  I 

K1=J1+I  5 

K2=J1+M 

C(Ki)=C( K1»-C(K2)^ST 
J1=J1+LL 
IJ  CONTINUE 
11  CONTINUE 
M1=M1+LL 
18  CONTINUE 
J1=C 

DO  <5  J=1  , LL 


1 F ( J-L«l J ) ) 1^,6,  14 

14  LKJ=LR (J ) 

J2= (LR J- 1 »»LL 
DC  13  1=1 ,LL 
K2-J2*l 
Ki=Jl+I 
S = C (K2  » 

c ( K^;»  =r.(  Ki ) i 

C(K1)=S 
13  continue 

LR(J)=Li‘  (LPJ) 

Lf  lLRJI  = LRJ 
IFU-LP(J))  14,8,14 
8 J1=JULL 
4 CuNTIMJt 
Rf  TURN 
END 
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VII.  EIGENVALUE  DETERMINATION 

The  subroutine  EIGEN(A,R,N,MV)  is  used  to  determine  whether 
lis 

G = Real(Y  ) is  a positive  definite  matrix.  This  subroutine  is 

•k 

taken  from  IBM  Scientific  Subroutine  Package  . It  can  only  be  used 
for  real  symmetric  matrices. 

MV  should  be  set  equal  to  1 (done  in  the  main  program)  for 
eigenvalues  only.  The  original  matrix  A must  be  stored  in  upper 
triangular  form  and  the  eigenvalues  are  outputted  on  the  diagonal  of 
A (therefore,  the  input  matrix,  A,  is  wlped-out  during  execution  of 
the  subroutine). 

Minimum  allocations  (for  eigenvalue  determination  only)  are 
given  by 

DIMENSION  A(N*N/2) ,R(1) 


IBM/360  Scientific  Subroutine  Package  (360A-CM-03X) , Version  III, 
Programmer's  Manual,  pages  164-165. 
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I 

ISUBRCUTiNE  EIGEN( AiR fNfMVl 

□ IMENSIGN  A(  14A) ,P(144) 

IF(^'V-l)  10  f 2 5 1 1 0 

10  IU=-N 

I DO  2C  J=1,N 

ig=iG+N 
DO  20  1=1 ,H 
I J=IQ-H 
R(  I J)=0.  0 
IF(I-J)  20,15,20 
15  R(  I J)=1.0 
20  CONTINUE 
25  AN0RR'  = 0.  0 
DO  35  1=1, N 
DO  35  J=I  ,N 
IF(I-J)  30,35, 3u 
30  IA=I*( J»J-J)/2 

ANGRM=ANaRM+A( IA)*A( lA)  ^ 

35  CONTINUE  I 

IF(ANORM)  165, 165, 

40  AN0RW=1.A14»S0RT(ANDPM) 

ANRf'X  = ANORM*'l  . 0E-6/FL0  AT  ( N ) 

IND  = 0 
THR=ANGF  M 

45  THR=THR/FLOAT(N) 

50  L=1 
55  M=L+1 
6U  MQ=(m*m-M)/2 
La=(L’t‘L-L  1/2 
LM=L+MQ 

62  IF(ABS(A(LM))-THR)  13.,b5,65 
65  IND=1 
LL=L«-LQ 

MM  = (i«  + |4Q 

X = 0.5*(A  (LD-AI-^M)  ) 

6a  Y = -A(LM)  / SORT  ( A(LM)«A(  LM  KXf'X) 

IF(X)  7^,75,7  5 
70  Y=-Y 

75  SINX  = Y/SQRT(2  .u«(  1.  SORT  tl  .0-Y>fY)  > ) ) 

SINX2=SI NX*SINX 
7a  COSX  = SQf-  T(1.0-SIMX2) 

CUSX2=CCSX-*C0SX 
SINCS  = SI  NX*CUSX 
ILG  = N*(L-1) 

IMG  = N’*(f^-l) 

DU  U5  I = 1,N 
I g={  I*-!-!  )/2 
IF(I-L)  8'J,115,aO 
30  IF(I-M)  35,115,90 
85  1^1=  I ♦MG 
GU  TO  95 
90  IM=MtIQ 

95  I F II -L)  ICO, 105, 105 
100  IL=I+LQ 
GO  TC  110 
105  IL=L^IQ 

no  x=A(m’*ccsx-A(i'")»siNx 
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A(  iy)=A(  IL)*SINX*-A(IM)*CUSX 
A(  lU  =X 

115  IF(MV-l)  12C  , 125,  l^J-j 
120  ILR=ILQ+I 
IMR=IMQ+ I 

X = R(  ILR)  *CaSX-R(  IMR)=»S1NX 
R(  IMR)=R  ( ILR)»SINX+k(IWR  l^'CCSX 
R( ILR) =X 
125  CCiNTINUF 

X=2.0*A(LM)*SINCS 
Y = A(LL)'‘CCSX2>-A(MM)»-SI  NX2-X 
X = A(LD*  SINX2*A(HM)»C0SX2^X 

A{LM)  = (A(LU-A(MN1)  ) « S I NC  A ( LM  ) ■*=  ( Cuo  X2-b  I 'J  X2  ) 

A( LL)=Y 

A(Mf»)=X 

130  IF(f'-N)  135,140,135 
135  M=M+1 

GO  TO  60 

140  IF(L-(N-D)  145,  150,  145 
145  L=L+l 

GO  TO  55 

150  IF(IND-l)  160,155,lo0 
155  INC=0 

GO  TO  50 

160  IF  (THR-ANRMX)  16  5,16  5,45 
165  IQ=-N 

DO  185  1=1, N 
I U=I Q+N 

LL=I>(  If  I-I )/2 
JQ=N*( 1-2) 

DO  135  J=I,N 
ja=jQ+N 

MM=J*( Jf J-J)/2 
IF  (A(LL) -A(MMI  ) 170,135,185 

17u  X=A(LL) 

A(  LL)  =A(  M)') 

A<N)^)=X 

IF(MV-l)  175,185,175 
175  DO  180  K=1,N 
ILR=IC+K 
IMR=JQ^K 
X=R(ILP) 

R( ILR) =R ( IMR) 

180  R( IKR)=X 
185  COMINUL 
RETURN 
END 


A 
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VIII.  MAIN  PROGRAM  AND  SAMPLE  INPUT-OUTPUT  DATA 

The  preceding  six  sections  have  described  all  of  the  necessary 
subroutines  (except  the  optimization  program)  to  generate  the  graphs 
contained  in  this  report  and  other  output  data  not  included.  The 
missing  optimization  subroutine  is  included  in  reference  [1].  This 
subroutine  is  a univariate  optimization  program  which  was  found  to 
be  very  efficient  (minimum  CPU  time)  for  optimizing  the  ratio  of 
quadratic  functions  normally  encountered  in  reactlvely  loaded  arrays 
(see  reference  [12]).  The  CALL  statement  for  this  optimization  sub- 
routine in  the  MAIN  program  (following  8 CONTINUE)  reads 

CALL  MAXU (N , NFP , N9 , STEP , BLOAD , PMER , YHS ) . 

All  of  the  calling  arguments  have  been  defined  previously  except  N9  and 
STEP.  N9  equals  the  number  of  complete  searches  (one  search  per  revised 
load  value)  and  STEP  is  related  to  the  size  of  the  Incremental  change  in 
the  load  during  the  searching  part  of  the  program.  Different  STEP  values 
should  be  tried.  Too  small  a step  will  require  more  searches  for  a given 
relative  maximum  while  too  large  a step  might  cause  a maximum  to  be  skip- 
ped past.  One  can  use  other  optimization  subroutines.  The  output  should 
yield  the  optimum  loads  for  maximizing  the  gain  corresponding  to  a given 
measurement  vector  P*"  (direction) . 

The  purpose  of  the  MAIN  program  is  to  obtain  input  data  and  call 
subroutines  to  generate  output  data.  Input  data  is  read  into  the  MAIN 
program  according  to 

READ  (1,201)  N,W,L,NEP,NA,N9,STEP 
201  FORMAT  (12,  2F5.2,  212,  113,  1E9.2) 

READ  (1,214)  (ANGLE(1),I=1,NA) 

214  FORMAT  (1013) 

READ  (1,216)  ER 
216  FORMAT  (1F7.2) 

READ  (1,216)  L3 
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There  are  N apertures  w-^ide  and  L-long  (L  Is  assumed  to  be  one- 
half  wavelength).  The  x-coordinate  of  the  T-th  aperture  is  X(I).  The 
feed  aperture  number  Is  NFP.  All  of  these  arguments  have  been  described 
previously  with  the  exception  of  NA, ANGLE, ER,  and  L3.  NA  is  the  number 
of  angles  for  which  the  gain  will  be  optimized.  This  integer  should 
correspond  to  the  number  of  angles  read  in  to  form  the  array  ANGLE(I). 

ER  is  the  relative  dielectric  constant  of  the  dielectric  filler  for 
the  waveguide.  This  number  can  and  probably  would  be  different  than 
the  mlcrostrlp  substrate  dielectric  proposed  for  the  reactive  loads. 

L3  is  the  length  of  the  waveguide  used  in  the  calculation  of  the  short 
positions  from  the  optimum  reactive  load  values. 

Minimum  allocations  are  given  by 

COMPLEX  PMH(N,37) ,PME(N,37) ,PMER(N) ,YHS(N*N) ,V(N) ,VCMR(N) 
DIMENSION  A(N*N/2) ,ANGLE(37) ,BL1(N) ,BLOAD(N) ,G1(37) ,G2(37) , 

G3(37) ,G4(37) ,NP(N) ,R(1) ,SHORT(N) ,VRMR(N) ,X(N) 

DO  loop  1 puts  kx  in  X.  DO  loop  2 calculates  the  eigenvalues 
lis 

of  Real  (Y  ) and  will  stop  the  program  if  any  of  these  values  are 
negative.  DO  loop  13  calculates  the  short  positions  (formula  19  in 
Part  I)  for  the  optimum  reactive  loads. 

The  sample  input  is  for  an  eight  element  aperture  array  with 
feed  element  aperture  number  4 and  aperture  spacing  0.29X.  The  gain  is 
optimized  in  the  ij)  = 50°  direction.  Loads  which  resonate  real  equivalent 
sources  were  used  for  initial  starting  points  for  the  optimization  pro- 
gram. (Note  loads  which  resonate  complex  equivalent  sources  could  have 
been  used  by  changing  the  CALL  BL0ADR(  ) statement  to  CALL  BL0ADC(  ).) 
The  MAXU  optimization  subroutine  went  through  six  complete  searches  for 
each  aperture  load  value.  Output  gain  values  are  given  in  5°  steps  start- 
ing at  0°  for  the  following  four  conditions. 

1)  Maximum  Gain  — Complex  equivalent  voltage  sources 

2)  Maximum  Gain  — Real  equivalent  voltage  sources 

3)  Maximum  Gain  — Using  reactive  loads  which  resonate  real 

equivalent  voltage  sources 

4)  Maximum  Gain  — Using  the  optimum  reactive  loads. 


c 
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COMPLtX  PMH(9,37),PME(9,37),PMFR(9»,YHS(40) 

COMPLEX  V(0I,VCMK(9) 

dimension  a ( I i^)  ,ANGLE(  37),HL1  ( 9»  ,!3LnAO(9  ) , G1  ( 37  ) ,G.^(  37  » ,G3  ( 37) 
dimension  G4(  37)  , NP(9)  ,P  ( 144)  , SHORT  ( 9 1 , VR>1  R ( 9)  , X ( 9 ) 
real  L,L2,L3 
INTEGER  ANGLE 

200  FORMAT (/ 3 X, 'N* , 6 X , • W • , 7X , ' L'  , AX , • WE P • , 9X , • NA • , oX , • N9 • ,7X, ' STEP' ) 

2 01  FORMAT! I 2 ,2E5.2,2I2,1I3,1E9.2) 

202  FORMAT  (/ 2 X,  1 2 , 3X , 1F5 . 2 , 3X,  IF  5 . 2, 3X , I 2 , 5X  ,I  2 , 5X , I 3 , AX  , IE 9 . 2 ) 

203  F0RMAT(9F5.2) 

20A  FORMAT (/ 3X, • X • , / / ( 3X  ,9F5 , 2 )) 

205  FaRMAT(/2X,'  Y— HALF  SPAC  E • / ( 3X , 5 E 1h  . 7 ) ) 

207  FCRMAT(/3X, 'APEPTURE  - -',I2,6X,'PM  =',2E1A.7> 

2C8  FORMATI//2X, • P - - MEASUREMENT  VECTOR  - - l-PLANE'/) 

210  FORMAT  (/ 3X, 'MAXI  MUM  GA  I N--CuMPU  X EQUIVALENT  VbL  T AGE  - -l  U J 1 V AL  E f.T  C 
lOMPLEX  VOLTAGE  FOR  MAXIMUM  GAIN  - -*,3X,13,'  UEokthS'/) 

211  FORMAT!/ IX, • MAXIMUM  GAIN  — REAL  EQUIVALENT  VOLTAGE  — cJ'JIVAlcN’ 
IREAL  VOLTAGE  FOR  MAXIMUM  GAIN  - -',3X,13,'  DtGKtLS'/I 

213  FORMATI/IX,'  LUADS--BL3ADS  - '13,'  UEGRFES'/) 

21A  FJRMATI10I3) 

215  FORMAT!/3X, 'PATTERN  AN  GL  E S • / / ( 3X  , i I 3 , ' LE35i.=i>'/n 

216  FORMAT! 1F7. 2) 

217  FORMAT  1/ 3X,  *P  ELAT  IVE  PERMITTIVITY  OF  UIELECTRIC  USc3  iN  wAVLGUli.r 
lELEMENTS  =',1F5.2) 

218  FORMAT !/ 3X,  »DI  STANCES  OF  SHORTS  F^R  rtAVEGOU'E.  LOAuS  - ',13,'  OEQC' 
lES'//» 

219  FORMAT  !// IX,  ' E-PLANE '//) 

220  FORMAT!/ 3X, 'GAI  N PAT'^’LRrjS  - -',3X,I3,'  UEGRr.ci'/) 

221  FORMAT  (/ 3X, 'MAXIMUM  GAIN  - COMPLEX  EQUIVALrNT  VuLT  AG  E ' / / ! 3 X , 7 F 7 . » 

U 

222  FORMAT  !/ 3X, 'MAXIMUM  GAIN  - REAL  EQUIVALENT  V OlT  AGE  ' / / ( 3 X , 7 F7 , 2 ) ) 

223  F0RMAT!/3X,'GAiri  - LPADcO  APERTUPES  - Ml  DA  L RcSuNANCc  L 0 AD  S ' / / T3  X , 

17F7.2  ) ) 

22A  FCRMATI/3X, 'GAIN  - LOADED  APERTURES  - OPTIMUM  LL  ADS  * / / ( 3X  , 7F  7 . 2 ) ) 
225  FORMAT !/ 3X, 'SHORT (',  1 I ,' ) = ',1F6.3,'  MM') 

220  F0RMAT(/3X, 'EIGENVALUE  ',12,'  = ',IE1A.7) 

227  FORMAT  !/ 3X, 'LENGTH  OF  APEPTURE  =',1F6.2,'  MM') 

WRITE!3,200) 

R£AD!l,2Ul)  N,W,L ,NFP,NA ,N9, STEP 
WRITE! 3, 2 02)  N,W,L,NFP,NA ,N9,STEP 
READ!1,21A)  ( ANGL E (I ) , I ^ 1 , NA ) 

WRITE!3,215)  ! ANGLE ( I ) , I = I ,N A ) 

REAC!1,216)  ER 
WRITE!3,217)  ER 
REAC!1,216)  L3 
WRITE!3,227)  L3 
MV  = 1 

PI=3. 191593 
PI 2 = 2.»P  I 
ETA=376. 730 
NN  = N*N 
Ni=N-l 

READ!  1,203)  !X( I ) ,I  = 1,N) 

WRITE!3,2U9)  (X!I),I=1,M) 

DO  1 I=1,N 
X! I)=X11 )*PI2 
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1 CQNTINUL 
L2=L 

CALL  YHbP(N,X,W,L2,YHS) 

*JSITE(3,2a5)  ( YHS(  J)  ,J=1  ,fiN» 

C POSITIVE  JcFINITt  TEST  FJf-  CCNOUCTANCE  MATi^IX 
JA=0 

Du  2 J=1 ,N 
JG=(J-L)*\ 

DG  3 I=l,J 

JA=JA+1 

JG=JG^l 

A( JA)=RL AL(YHSIJG)  ) 

3 CGNTINUf 

2 CCNTINUF 

CALL  EIGEN(A,R ,N,MV) 

J2  = l 

DJ  4 I=1,N 
WRITE(3,226)  I,A(J2) 

IF(A(J2) .LT.O. ) GC  TL  5 
J2=J2*I>1 

4 CONT  INUE 

CALL  PATHSPINtX.Lf W,PME,PMH» 

DO  6 Ii=lfMA 
I = ( (ANGLE  II  1)*36)/1«C+1) 

J1=ANGLE  (ID 
KRITE(3,2C8> 

DJ  7 J=1 ,N 
PMER( J»=PPE(JtI) 

WRITE(3,2C7»  J,PMEP(J) 

7 CONTINUE 
WRITE(3, 219) 
wmTF(3,210)  J1 

CALL  MAXGCV(N,PVEP ,YHS ,VCMR) 

WRITE(3,211)  J1 

CALL  MAXGRV(N f PMERiYHS fVRMR) 

WRITE(3,213)  J1 

CALL  BLOAURIN ,NFP,BLCADtVRMR , YHS ) 

DO  8 J=1  ,N 
BLK  J)=ELuAD(  J) 

8 CONTINUE 

CALL  MAXU(N»NFP,N9,STEP, BLOAU,PMER,YHSI 

DO  9 K=1  , N 

V(K)  = ( 1. ,0.)*VR^R(K) 

9 CONTINUE 
WRITE(3,  220)  J1 
DO  10  M= 1 ,37 

DO  11  )=1,N 
PMER( J)^PPE( J,P) 

11  CONTINUE 

CALL  FUNC  TA(N,GA1 N,PMER, VCMP ,YHS) 

G1 IP)=GA IN 

CALL  FUEjCTA(N,GAIN,PMER,V,YhS) 

G2(P)=GA IN 

CALL  FUNCTE'(M  ,NFP  ,3L  1 , GA  I N , P ME  P , VMS  ) 

G3(M)=GA IN 

CALL  FUr'CTH(N,NFP,BLOAD,GAlN,PMEK,YHS) 

G4(M)=GA IN 
10  Continue 
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WRITE(3,221)  (GUM)  ,M=1, 37) 

WRITE(3,222)  ( G2 ( M ) , M= I , 37 ) 

WRITE(  3,  223)  ( G3  ( M ) , M=  1 , 3 7 ) 

WRITE(3,22<»)  (G^(M)  ,M=l,37) 

12  = 1 

OQ  12  I=ltNl 
IFU.EQ.NFP)  I2=I2+1 
NP(  I )=I2 
12=12+1 

12  C3NT1NUE 
WRITE(3,21d)  J1 
DO  13  J=1,N1 
J2=NP( J) 

SHORT ( J2 ) =ATAN (SQRT(ER-l.)/(-BLOAO(J2)+ETA))/((Pl*SUkT(tR-l.))/l  ') 
WRITE(3»225)  J2fSHGkT(J2l 

13  CONTINUE 
6 CQNTINUE 
5 CONTINUE 

STOP 

END 


$DATA 


N 

L 

NF  P 

NA 

N9 

STEP 

8 

0.0  1 

0.  5C 

4 

1 

49 

0. 10 E -06 

PATTERN  ANGLES 
50  DEGREES 

RELATIVE  PERMITTIVITY  CF  DIELECTRIC  USED  IN  WAVEGUIDt  ELEMENTS  = i.-<P 
LENGTH  UF  APERTURE  = 51.72  MM 


X 


-l.l6-u.87-J.58-C. 29  J.OC  0.29  J.58  ".87 


Y — HALF  SPACE 

0.41 1908 3E-04  0 . 2343 4 55E-04  0 . 17 80 1 1 8E -U 4-u . 1 8b 5 384t -3 4-C . 12  3321  It- u4 
-0.1068729E-L4-0.6214996E-Q5  0 ,9807  1 1 3E -!•  5 . tt32d  394F-J  5 0 . 324461  7E-.  d 

0.1 13210  7E-C5-0.  7153  6 54E-05-0.606448  5r-05  0 .4 1 98 1 54E -J 6 C,  15 57 1 l8b-E 5 
0.4996528E- J5  0 . 1 78l  1 18E-04-C  . lb  85  384E 4 u . 4 1 1 9(.  8 3E -D  4 0. 234j455l-^h 
0.17801 14E-04-0. 1 885 384E-04-0 . 12 3321 3 ^ -04-0. 10 68 726E-D 4-0 . 62 149b4E-u5 
0.9807  145E-05  0 .8328  398E-05  ...  32  4461 7E 5 ^ . 1 1 32  I 4F -3  5-^ . 7 1 5i634L- -5 
-0.6064482E-05  0. 4 1984 76E-06-0 , 12 3321  1 E -U4-0 . 1 068 729E-J 4 0 . 1 7oj 1 14£- J4 
-0.1885384E-04  0.41  19083E-04  0 . 2343455E -o4  ^ . 1 78t  1 18F-3 4-u. io853b4L-u4 
-0.123321  1E-04-0.1066  729E-04-0.6214996F-0  5 0.980  7 1 13E -3  5 0.  83  28  3'5' 4E-^  3 
C.3244617E-G5  0 . 1 1 32  104E-05-0 . 71  53654E-.  5->'  .fc^l4996E-35  v8u7  1 1 3L- . 3 
-0. 1233 21 3E-04-0. 1068  7 26E-04  0.1  7 80 1 1 8E-C4-0 . 1 885  384E-3 4 0.41 l9u6  3E-u4 
0.234  3455E-04  0. 17dt  1 18E-04-0 , 18  85  384F 4-'' . 1 2 33 2 1 It-U 4- < . lo 6o  7t 9t- ^ 
-0.6214996E-05  C .9807 1 13E-05  0. 8328398E-05  U. 3244617E -3  5 0. 832a394c--5 
0.32446  17E-05-0.6214  9 84E-0  5 0 .9807  14  5E-i.  5-0.  12  332 1 IF.-J  4- 0 . 1u667«:9l-s.-* 
0.1  7801 18E-04-0.  1885  384E-04  0 . 41  190a3r -U4  U.  ? 343  455E-J  4 0 . 1 78- 1 1 o t-^ -» 
-C.  1885  384  E-C.4-  J . 12  33  2 I lE-04-:.'  . lu68  729F  -t  4-' • . 62 1 4 9d4E  -3  5 o.9ou7  l4  5r.-^  5 
0.113210  7E-05-0. 7153654E-05  0 . 8328398E -0 5 C. 32446l7t -3 5-G. 62 14996L-U 5 
C.98C  71 13E-05-0.  1233 2 1 lE-G  4-0 . 10  68  729F-r  4 . 1 78v  118E-j4-vj.  Ibd53b4c-.i4 
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0.411VC8  36-C4  0.2i43^55E-04  0 . 1 78U1 1 8E-04-C. 1 6o5  3«4E -3 4-0. l2  33213c-04 
-0.106  8 7^6E-(  4-0  .6064  4 85t'-05  .)  .4198  154E-.  6 u.  1 132  K'4E-3  5- v . 7 1 53654c- -/ 5 
0.8  328  39  4E-05  0.  32  4*t  6 17E-0  5-0 . 62  14  996E -0  5 0 .9  c»07  1 1 3E -D  5-0.  12  3321 1 c-U*t 
-0.1068729E-04  0.178*  1 18E-04-C' .1885384P-:  4 .411Vj83L-34  0.  ^3434551  - .'4 

O.l  78CI  14E-04-0.  1885  3 84E-0  4 0 . I 5 57  1 1 8E-vi  5 r . ^.996  528E -D  5-f;.  60644<}2E- J 5 
0.4198476E-06  0 . 1 1 32  1C4E-0 5-C  . 7 1 53654E -0  5 f . 8 32o 398E - J 5 *j.  J^446 1 7E- - 6 
-0.6214984E-C5  C .9  807  1 4 3E-Q5-0 . 12  33  2 1 3E -0  4-0 . 1(j60  726F -3  4 0 . I7oo  1 l4  E- .»4 
-0.1885384F-04  0 .4 1 19083E-04  0 . 234345 5F -04 


EIGENVALUE  1 = 
EIGENVALUE  2 = 
EIGENVALUE  3 = 
EIGENVALUE  4 = 
EIGENVALUE  5 = 
EIGENVALUE  6 = 
E IGENVALUE  7 = 
EIGENVALUE  8 = 


0. 7963572E-C4 
0. 748 1415E-C4 
0.59 J2007E-C4 
:.55 14j34E-C4 
0.4829286E-04 
0.118  1671E-04 
0.  7b98826E-i'6 
3.1440352E-l7 


P - - measurement  VELTCR  - - E-PLANE 


APERTURE 

— 

- 1 

PM 

= 

0.21 9439  76 

-w  1- 

V • 

7996443E 

- 'j 

APERTURE 

- 

- 2 

p y 

= 

0 . 745  1952F 

Ov*  — 

0. 

29  /c5  52E 

APERTURE 

- 

- 3 

PM 

= 

0. 5578344E 

jO 

0. 

5 ^3526E 

*..0 

APER  TURE 

- 

- 4 

PM 

-0. 311 1879L 

OL 

0. 

7369359c 

V^O 

APERTURE 

- 

- 5 

PM 

= 

-0. 7999454E 

Jo 

0. 

C 'oc  C ^ OL 

LC 

APERTURE 

- 

- 6 

py 

= 

-0.31  118  79E 

^0- 

0. 

7369359E 

cO 

APERTURE 

- 

- 7 

PM 

0.5578344E 

uO- 

0. 

5733526c 

-.0 

APERTUP  E 

_ 

- 8 

PM 

— 

0. 745  L952E 

03 

J. 

29.J6552E 

>jO 

E-PLANE 


COMPLEX  gain  =10.24 

COMPLEX  ECUIVALEM  Vl.  L TAGE  FGF.  MAXIMUM  CLmpLcX  GAIN 

-0.5144u69E  05  0.6518969E  05  0.129431aE  06-0.  lobS^lbL  3 6-0.  1 7.:352.  L 
0.19353C8E  Ci6  C.1972  7 70E  06-u  . 28  37  579r  v 6-u . 1 o45980F  36  >,'.^9ti6olL  .6 
0.llll99oE  C6-0.23411C3E  06-O.4962431f  05  0.1blH'»*«9t  J6  0.3u8.;o3aE  u'~> 
-0.6dfc73SA£  C5 


MAXIMUM  PtAL  UAIK  = 5.75 

PEAL  EQUIVALENT  VCLTAGE  FOP  MAXIMUM  PEAL  GAIN 
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~0.1347424E  L6  u.293254dE  06-^.4644^18^  i'6  ( .62531p9£ 
G.4645o5tiE  06-C . 2933  5S6E  06  0.1347722K  u6 


j6-U.o253364L 


j o 


LOADS — BLLADS  - 50  DtGPEES 


bljad-pescnating  keal  sources 


-0. 1774331E-J4-0.3C  73  0 09E-&4-0  .37  55  36  3E-C4  C.o-.  JCOc  UF  OO-C.  329 j 3 1- . 4 
-0.375567«E-04-0.3073015E-04-u. 1774960t-04 


HI 

= 0.140 99 69E 

C3  0. 

JCCOODOE 

00 

H2 

= -(-.2527136E 

<6-0. 

l9u47l5E 

07 

H3 

= -0. 25272  lie 

C6  0. 

1904705E 

07 

H4 

= 0.1545716E 

12-0. 

3 3952c(.'E 

05 

H5 

= 

t’.111956dE 

C9 

H6 

r 

-0. 1417972E 

C5 

H7 

= 

O.2964023E 

C3 

H8 

= 

0. 3423529E 

C6 

H9 

-0. 11 70265E 

(,5 

U1 

-L . 1 1 1 6b07E 

(..5  j.2477C63E  04 

DEN  = 0. 1349386E  03 
GAIN  = 2.32 


MINIMIZATION  OF  Y = F(X)  WITH  RESPECT  TO  HH  LOAD 

X(l)  = 0.5635927E  l5 
Y(  1)  = 0.*i31e356F  00 
S2  = C.  1769241E  l5 


PRELIMINARY 

SEARCH  ON  the 

NEG 

ATIVE  SICE 

0.5635927E 

05 

0.5635927E 

U5 

n .56359<i7E 

0 5 ..56  3 592  21- 

0 5 

-..5635913L 

5 

0.5635899E 

05 

O.  5635071E 

05 

0. 5635816E 

J5  <j.5bj57juE 

05 

0.  5o354S9t; 

^5 

0.5635012E 

05 

0.43l635e£ 

00 

C.4316353E 

00 

0.4316353E 

vO  0.4  3lb  35  7E 

OO' 

0. ^3 lo3o3L 

C.4316371E 

00 

0.4316385E 

CO 

0.43 1641 7E 

:0  ^'.43i6483t 

O'; 

0.  43  160  1 £. E 

0.4316873E 

uo 

PRELIMINARY 

SEA 

,RCH  ON  THE 

POSITIVE  SIDE 
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0.5635927E 

05 

0. 56  35  0 12E 

05 

0.56 35 93  IE 

u5 

0.5635931E 

D5 

0,5635936l 

5 

C.5635945E 

05 

0.5635959E 

05 

0. 5o35987F 

C'5 

j . 5o  36  o42E 

0 5 

0.  5b  3b 1 58  c 

0 5 

0.5636389E 

05 

0.5636847E 

05 

U.4316356E 

-)0 

C.4  316  8 7 3E 

00 

J.43  1635:>t 

'■.4^16  353!: 

O'. 

6.931634  6;: 

y ./ 

0.4316343E 

00 

C.4316335E 

00 

0.431632.^6 

5 1 1. 2 896 

3 ) 

w. 4j i62t  2 L 

0.4316092E 

00 

C.4315833E 

00 

FINAL  SEARCH 

0.56359276 
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XI  1)  = 0.325-^1  J9E  C5 
Yin  = 0.377a77lf  .lO 
SZ  = 0.  I4l2<*73t  C5 


MINIMIZATION  OF  V = F(X)  ^ITH  FESPECT  TO  3'^H  LOAu 

Xm  = 0,266285eE  15 
Y(l»  = 0.26568»7E  i)0 
SZ  = 0. I Ae587tiF  C5 


MINIMIZATION  OF  Y = F(X)  WITH  FlSPcC'’  TF  5TH  LCAU 
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MINIMIZATION  OF  Y=F(X»  WITH  FESPECT  TO  oTH  LuAU 
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D1 STANCES 
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